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Abstract 

The sphaleron rate is defined as the diffusion constant for topological number N^g = J ^^^f . 
It establishes the rate of equilibration of axial light quark number in QCD and is of interest both 
in electroweak baryogenesis and possibly in heavy ion collisions. We calculate the weak-coupling 
behavior of the SU(3) sphaleron rate, as well as making the most sensible extrapolation towards 
intermediate coupling which we can. We also study the behavior of the sphaleron rate at weak 
coupling at large Nc- 



I. INTRODUCTION 



It has long been known [l| that Yang-Mills theory possesses multiple classical vacua, 
distinguished by a topological number, the Chern-Simons number. No continuous sequence 
of infinitesimal gauge transformations can bring one of these vacua into another; instead a 
"large" gauge transformation of nontrivial topology is needed. The topological character of 
the gauge transformations is ensured by the nontrivial homotopy of the gauge group, e.g. 
TTsiG) = ZioiG = SU(iVc), > 2. 

This topological nature of the vacuum has important physical consequences because of the 
chiral ( AB J) anomaly j2|, [sj . The would-be conserved current associated with a left-handed 
fermion in the fundamental representation is^ 
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Within the Standard Model, this has a few very interesting consequences. The SUl(2) weak- 
hypercharge group couples only to left-handed fields, so the anomaly gives rise to baryon 
and lepton number violation While under ordinary conditions these have essentially 
no physical consequences, at very high temperatures where electroweak symmetry is "re- 
stored," baryon and lepton number violation becomes efficient jsj-Ql, opening the possibility 
of electroweak baryogenesis The SUc(3) color group couples to equal numbers of left 
and right handed fields^ and so no global symmetry violation is associated with it. However, 
the up and down quark fiavors are extremely light, so the anomaly is the most efficient way 
of interchanging left with right handed first generation quarks. 

The equilibration of right- and left-handed quark number is irnportant to baryogenesis, 
and therefore the strong sphaleron rate is too; see for instance |9-12|. The same applies 
to heavy ion collisions, which also appear to involve the production of a hot quark-gluon 
plasma 13|. Indeed, recently Kharzeev, McLerran and Warringa have argued that topolog- 



ical fiuctuations may play an important role in creating local fiuctuations or imbalances in 
left-handed versus right-handed quark number IJ, |l5|, which could have very interesting 
experimental signatures 16 . 



The relation between real (Minkowski) time topology change and Euclidean topological 
susceptibility is subtle (see [l7| and the discussion below). Therefore the Minkowski rate 
of topology fiuctuation cannot be easily inferred from Euclidean quantities and needs to 
be separately studied. There is a vast literature on doing so for the SUl(2) sector of the 
Standard Model 

m 

■ lisl-lisl. and it would be fair to say that the efficiency of topology change 
is well understood for this group and at the relevant coupling. But SUc(3) has been much 
less carefully studied. At this time, there are two quite different parametric estimates [ol, [2^ 
and one rather crude calculation 3^ which was made before our modern understanding of 



the parametric behavior was developed 2J, |25 



recent resurgence of interest in topology change in the strong interactions [IJ, [15 



Therefore, particularly in the light of the 

we think 



^ Our conventions are that g^^, = 77^^ = Diag[ — + and that capital letters X,P are 4- vectors with 
t = x*', oj = and x,p are the 3-vector components. Lower case x,p are reserved to denote |a;|, 

^ or equivalently, to equal numbers of left-handed fields in the fundamental and antifundamental represen- 
tations 
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it is timely to revisit the question of the "strong sphaleron rate," the (equihbrium, thermal) 
rate of topology change in the group SU(3). As we will discuss momentarily, we currently 
only know of tools to do so at weak coupling, so we will be forced to restrict our attention 
to weak coupling. However we will also make the most reasonable estimate we can of the 
rate at intermediate coupling. 

For the impatient reader, we summarize our results here: the rate of topology change per 
spacetime volume is defined as 

Tsphai = hm (^cs(^) -A^cs(O))^ ^ J / J^F;,Fr(X) ^,F:,FfiO)) (1.2) 



i— >oo 



where the t — t- oo (zero frequency) limit is to be taken with some care, see below. In 
the limit of high temperature (and hence of weak coupling), the temperature and coupling 
dependence for SU(3) is (as = fi'^/47r) 

Tsphai = (132 ± 4) (^^^ (^In ^ + 3.0410^ a^T' , (1.3) 

7= In ^ + 3.041 V (1.4) 



47r \ 7 



K = g g'T' = —f^g'T' . (1.5) 

Note that 7, the mean rate of color randomization, is defined self-referentially, and should 



be found self-consistently [28 



The A^c scaling is surprisingly simple; investigating Ac = 2, 3, 5, 8 we find that 

r,phai = (0.21 ± 0.01) (^^) (in ^ + 3.0410) ^^^{KafT' (1.6) 
with no further A^. corrections visible within our error bars. The functional form and even 



the numerical value are in good agreement with the conjecture of [29 . 

The formal range of validity of the above results is very narrow; not only the coupling ctg, 
but even the inverse log of the coupling l/ln(l/Q;s) must be small. At larger couphng we 
can estimate the sphaleron rate based on the diffusion rate of a lattice-regulated, classical 
field theory with the same magnetic-field damping rate as the leading-order result we find 
in the quantum theory. This approach is known to be correct in the weak coupling limit 



2J, |25|, |31| but at larger coupling values it amounts to a crude estimate; the dynamics cease 
to be well described by classical field theory, and the classical lattice model sees more and 
more "lattice-y" dynamics artifacts. We can also only push the lattice to a certain level 
of coarseness before it becomes impossible to establish any topological character for the 
fields. As shown in Table UTTt this allows access to as ^ 0.1. The leading-log results work 
reasonably well for as ^ 0.03 but break down after this. And the classical framework starts 
to break down above about as = 0.1. 

The remainder of the paper is organized as follows. In Section [TTl we review the definition 
and utility of the sphaleron rate, its relation to Euclidean correlation functions, and the 
extreme difficulty of its determination from lattice Euclidean data. Then Section UTTl reviews 
the effective field theories valid at weak coupling and presents the details of our numerical 
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implementation and calculation. Section IIVI presents our numerical results for the weak- 
coupling sphaleron rate and Section |V] presents the best extrapolation towards intermediate 
coupling we can achieve. We discuss applications in Section IVII 



II. SPHALERON RATE: GENERALITIES 

Here we review the definition and physical relevance of the sphaleron rate, and its relation 
to Euclidean correlation functions. Nothing in this section is new, it is included as a refresher 
for readers familiar with the issues and an introduction for readers who are not. 



A. Definition of the sphaleron rate 

As mentioned above, the axial current associated with a quark species, 

Jl,, = qi'l^q (2.1) 

is not conserved. There are actually two contributions to its nonconservation: the quark 
mass, which explicitly breaks the axial symmetry, and the chiral anomaly:^ 

a, J^,^ = 2m,g75g - 2 J^F^.^r (2.2) 

(where the factor of —2 difference from Eq. f ll.ip is because is the right-handed minus 
the left-handed currents). In real- world QCD the up and down quark masses are very small. 



m„ rf ~ 2.0,4.5 MeV j32| and so in many situations one may neglect the explicit current 
nonconservat ion . 

Neglecting the quark masses, the amplitude to evolve from state \ipi{ti)) to state \'ip2ii2)) 
times the change in axial number is (in the Heisenberg picture) 

dt J J^F;,Fr(a;,t)|^i) . (2.3) 

The probability of the process times the square of the change is 

p{AQ^,,r = 4 1 d'xd'Y{i;^^^F;^Frmi^2)m^^ ■ m 

Summing over final states, \'ip2){'^2\ = Ij we find 



The mean squared change for a general density matrix p is found by replacing \ . . 
with (Tr p...). Note that the operators here are not time-ordered; it is the Wightman 
correlator 



^2 „2 



G>^(X, Y) ^ { 4-,F;^Fr{X) 4-,F:,Ff{Y) ) (2.6) 



From now on wc specialize to vcctorlike theories with fundamental representation matter. 
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which is relevant here. 

Eq. (12. 5 P holds for arbitrary density matrices including nonequilibrium density matrices. 
Here we will more modestly try to evaluate it for the thermal ensemble. Since the thermal 
ensemble is 4-translation invariant, the mean-squared QA,q change is extensive in 4- volume 
and it is more useful to define a rate of change per 4-volume: 

((A(5A,(j)^)thcrm. = 4V^t Fsphal , 



sphal 



j d'XG>^{X,0). (2.7) 



This is the same as the zero frequency and momentum limit of the momentum-space Wight- 
man function 

G>^{P) ^ I d^Xe-^^-^G>^(X, 0) . (2.8) 

In defining the sphaleron rate in this way we have tacitly assumed that the generated 
{QA,q) value can persist without influencing the subsequent evolution of the system. Of 
course this is not true; a net abundance of QA,q leads to chemical potentials for left and 
right handed quark number which will then energetically bias future topology change to 
allow that quark number to relax. Defining Qr = X]g Qm the total axial light quark charge, 
standard fluctuation-dissipation arguments [oi [s^ show that 

^ = -g^MZlEEi^ ^ _Q ^^ (2.9) 

where xq is the susceptibility for and in the last line we used the leading-order perturba- 
tive estimate xq = N{ N^T"^ /3. What this expression means is that any fluctuation-induced 
nonvanishing value for will bias all subsequent topology changing processes, causing 
to relax back to zero with fluctuations set by the susceptibility. This means that a fluctu- 
ation in topology of one sign will eventually be balanced by a fluctuation of the other sign 
with a time constant r ~ 2TxQ/(4A'^f ^Fsphai)- To define the sphaleron rate we implicitly 
assume that this time scale is long compared to any microscopic time scales, and we implic- 
itly cut off the time integration in Eq. (12. 7p on a timescale shorter than r but longer than 
microscopic time scales. At weak coupling the longest microscopic time scale is ~ l/«s^ 
while the relaxation time scale defined above is ~ l/^gT (as we shall see), so the separation 
of scales exists. The separation also exists at strong coupling in large Nc A^=4 SYM theory. 



where the sphaleron rate is Egphai = {g'^N^YT^ /256ti^ ~ [3J], while xq ~ N^T"^. There- 
fore T ~ /T is large in the large number of colors limit. ^ However for Nc = 3 ordinary 
QCD at T ~ 200MeV it is by no means clear that such a separation of scales exists, and it 
is possible that one cannot really define a sphaleron rate per se. 



B. Relation to Euclidean correlation functions 



We saw above that the sphaleron rate is set by the zero frequency and momentum limit of 
the FF Wightman function. Naively one might think it is easy to determine the sphaleron 



^ The calculation is only valid in the limit A^c ^ {q^^c) ^ 1- In this limit the relaxation time is paramet- 
rically large, while all microscopic timeseales are ~ T. 
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FIG. 1: Integration contours in complex t plane for and (the latter for 

lm.{u) = 27rn only). Contour deformation turns G^, but not G^, into G^. 



rate from Euclidean correlation functions, which can be measured on the lattice. Indeed, in 
the late 1980's there were some misconceptions that such things should be possible, which 
were quashed by a paper by Arnold and McLerran 17|] . Here we remind the reader how the 
sphaleron rate is related to Euclidean correlation functions, and why it would not be easy 
to determine it from Euclidean lattice data. Though the remarks here are rather generic to 
two-point functions and are well known, we review them in a little detail because we are not 
aware that they have been clearly discussed in this particular context. 
The Wightman correlation function defined above is 



G>^(X)^Tre-^^e^^*J^F;,Fri 



while the Euclidean correlation function is 



Gl,iT,x) ^ Tre-^^e^^;^F;,Fr(a.)e--^ 



327r 



327r 



(2.10) 



(2.11) 



defined on the range < r < /3 and satisfying periodic boundary conditions, G^{(3,x) = 
G^{0, —x) = G^{0, x).^ It is clear from the definitions that G^{t, x) is the analytic contin- 
uation T ^ it of G^{t, x). In particular the equal-time value of G^ equals the equal r value 
of G^ 

But what we want is the zero Minkowski frequency limit of G^ , and the frequency trans- 
forms of the two functions are not the same. Instead, the frequency-space transform of the 
Euclidean function continues to the retarded correlation function 



G ppit, x) 



m {GlJt,x)-G<Jt,x)) 



■ FF^' ' FF 

The reason why is illustrated in Fig. [1] first note that the frequency domain form of G^ 



(2.12) 



G-iuj) 



dte"^'(G>{t)-G<{t)) 



(2.13) 



is trivial to continue to complex u with Im w > 0. For purely imaginary uo = iu^ with 
f3ujE = 27m it is 

/•OO /"OO 

G^(^i^^)= ^^^i^BitQ>^^-^ _ ^^^iu^itQK^^y ^2.14) 

Jo Jo 



Note that FF is a pseudoscalar. G{—x) can be related to G{x) by a rotation; parity is not involved. 
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But the KMS condition is that 



= Tre-^^Fe-^^e'^^e^^^FFe-*^* (inserting e±^^) 

= Tr e-^^ e^^e'^'FFe-'^'e-^^FF (cychcity of the trace) 

= Tr e'l^^ ^i{-iP+t)H p p^-i{-i(3+t)H p p 

= G>^{t-zP). (2.15) 
Substituting into Eq. fl2.14p and noting that e'"^'^^ = 1, we find 

POO poo—ifi 

Ciico^) = / rfte*'^E**G>(t) - / dte'^^''G>it) 

Jo J-il3 

= / dte'^^''G>{t) = -i / rfre*^E"G^(r) = -iG^M (2.16) 
Jo Jo 

the Fourier (series) transform of the Euchdean function. 

The conclusion is that, while Euclidean techniques determine G^{t = 0, x) in a straight- 
forward way, the zero Euclidean frequency correlation function [the Euclidean topological 
susceptibility] corresponds to the zero Minkowski frequency, retarded function, not the zero 
frequency Wightman function. That is, the Euclidean correlation function continues to a 
Minkowski function with the wrong operator ordering. 

The Euclidean function does have some relation to the Wightman function, through the 
Kramers- Kronig relations. Namely, defining 

a^pit) = G>{t) - G<{t) a^^iu) = (1 - e~^^)G>{uj) ~ PujG>{u) , (2.17) 

the retarded function at arbitrary complex frequency u is determined by cr{uj') at real fre- 
quency by 

G-Ju)= r^^Fii^^ ^2.18) 
Evaluating for imaginary (Euclidean) uj and inverse transforming to r, 

doj' ^FF cosh [uj\t - /3/2) ) 

Unfortunately, what we can measure on the lattice, G'^^(r), determines an integral over 
what we want, app{(jj'). To proceed further we would need some method to invert this 
integral relation. This problem has frustrated the evaluation of similar equilibration rates 



such as the electrical conductivity and the shear viscosity [35 



We now argue that the problems encountered in continuing FF from lattice data are 
more severe than for continuing the correlation functions needed for conductivity or viscos- 
ity. The problem is that there is no local^ lattice definition of the topological charge density 



^ By local here we mean that the operator has compact support which is independent of the lattice con- 
figuration, i.e. it stretches for a fixed finite number of lattice spacings. It may be possible to meet our 
requirements with an operator which falls off exponentially with lattice distance, but only if the exponen- 
tial falloff rate is conditional on the smoothness of the lattice configuration. 
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operator -^:;^F°;u^b'^ which satisfies the property that J d^X^^F^j^Fj^^ over a compact 4 



volume is topological (for a discussion see 36[). There are definitions of / df^X^^F^^F^ 



which are topological (given some mild constraints on the smoothness of lattice fields) (37| . 
That means that there is no problem in principle with evaluating G^^p,{(jj^ = 0,p = 0) or 
[given real-time lattice configurations] G^p,{u = 0,p = 0). But any attempt to compute 
J d^x G^p{T, x) relies on measuring FF on a 3-D rather than 4-D surface, and so will neces- 
sarily be contaminated by either the non-topological or the non-local nature of the operators 
used in the evaluation. We believe that this additional complication will make Fsphai even 
harder than other transport coefficients to extract from Euclidean lattice calculations. 



III. EFFECTIVE THEORIES AND COMPUTATION 



As reviewed in the last section, the sphaleron rate is an intrinsically Minkowski quantity 
which cannot be evaluated with Euclidean methods. Because iSr^F^,F!^'^ is a total deriva- 
five, the sphaleron rate is a strictly nonperturbative quantity which vanishes to all orders 
of perturbation theory. Therefore its evaluation requires nonperturbative, Minkowski tools. 
These are in very short supply. To our knowledge, we have only two approaches to measuring 
the sphaleron rate: 



In certain gauge theories which possess gravity duals (generally highly supersymmetric 
theories) there are holographic methods which allow the sphaleron rate to be related 
to string theory correlators. These can be evaluated on the string theory side in the 
strong-coupling, large A^c limit of the gauge theory, see 3J] . The results are suggestive 
but do not apply to ordinary QCD. 

In the weak coupling (high temperature) domain, the only nonperturbative physics 
which is not exponentially suppressed is essentially classical Geld physics. Therefore it 
should be possible to evaluate the sphaleron rate by studying classical field dynamics 



as an effective IR description [18|, [38 



In this work we will follow this second approach. 



to its logical conclusions for the electroweak SUl(2) sector [20H26 
what that literature has found. 



This method has already been pushed 
Here we briefiy review 



A. Bodeker's Effective Theory 

Naively, the infrared fields in a weakly coupled quantum field theory of massless degrees 
of freedom should obey the classical equations of motion, in our case 

with the initial field values drawn from a statistical distribution of classical fields weighted by 
the Boltzmann factor exp{—E/T). Bodeker has shown [3^ that the first corrections to this 
picture enter at order h"^ except for the possibility of large (loop) corrections in the ultraviolet. 



8 



That is, is a composite operator so an infrared field component could arise from the over- 
lap of UV fluctuations with almost the same wave number, e.g. J^{k) = il}{—p)'^^ip{p+k). 
These fields can in turn be correlated with past and present IR gauge fields, leading to ef- 
fective modifications of the infrared gauge dynamics. Such large modifications of the IR 
dynamics do indeed occur; they are just the Hard Thermal Loops of Braaten and Pisarski 



40, 41 



The sphaleron rate depends in particular on the behavior of transverse (magnetic) fields 
in the deep infrared. In this case the Hard Thermal Loops have a particularly simple 
interpretation; the spatial current in Eq. (13.11) is determined by Ohm's law and the color 



conductivity; writing noncovariantly, the spatial components read [25 



D F ■ 



D,E, 



-Ji 



-cjE, + 



(3.2) 



with a the color conductivity and Q a Langevin noise source which maintains thermal 
equilibrium. In general the color conductivity is frequency and momentum dependent, but 
for the most infrared fields of relevance to the sphaleron rate, "color randomization" means 



that a should be essentially a constant of value |2J, |25|, [28 



a 



1 



4tt 



In + 3.0410 



2Nr 



6 



7 

—g^T^ . 



0{l\ 



(3.3) 
(3.4) 



(Physically, the conductivity a = m'^/S'-f is a polarizability times a mean free path. The 
factor of m'^/S is the polarizability and I/7 is the mean free path for color randomization.) 
This conductivity is parametrically large and allows one to neglect the DfEi term in Eq. (13. 2p 
relative to the aEi term, yielding the effective description 



DF- 



-aEi + Ci 



dtCiix, t)Cj{y, 0) = 2aTS,,S%x - y) 



(3.5) 



The normalization of Q ensures that the fields are distributed according to the classical 



thermal ensemble, as required by the fluctuation dissipation theorem 25 



Note that in finding this effective description and in particular that a is a constant, 
we have made an expansion (to second order [2^) in hi{l/g) the logarithm of the inverse 
coupling. While the numerical coefficient 3.0410 stabilizes this expansion somewhat, it 
is a rather marginal expansion and one should expect rather large corrections at physically 
interesting couplings. Nevertheless we first pursue this approach because the effective theory 
involved, Eq. (13. 5p . is a Langevin equation with very good mathematical properties. In 
particular, unlike Eq. (13. ip . Eq. (13. 5p possesses dynamics which have a well defined limit as 



the UV regulator is taken to infinity 25|, |27|. Therefore it can be implemented on a lattice, 
and the lattice spacing a — > limit may be taken to obtain a well defined limit for Lgphai in 
this effective theory. 
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B. Numerical implementation 



Our implementation of SU(3) gauge theory on a lattice is standard, and we measure 



Chern-Simons number using the algorithm of 421], which extends easily to general SU(A'"( 



We do not follow the implementation of Langevin dynamics from [27j, however, so we will 
make some remarks about how we do implement Langevin dynamics. First, it is convenient 
to rescale the time variable in Eq. fl3.5p to conventional Langevin time 

Tl = t/a (3.6) 

which, note, has units of time squared. In terms of this Langevin time the fields evolve 
according to conventional Langevin equations; in Aq = gauge, 

f)]-f f 

d^A, = -D,F^, + C = + C, J dT^Cir,, y)C,ix, 0) = 2T6,j6'ix - y) . (3.7) 

The sphaleron rate Fl = / c/tl / d^xFF^Ti^, x)FF{0, 0) obtained in this theory is related to 
that obtained using Eq. (13. 5p by a factor of a. 

Any numerical implementation of Eq. (13. 7p requires discretization of Tl into finite steps. 
The challenge is that larger steps result in a faster algorithm but may have larger numerical 
errors; in particular the sphaleron rate scales as T"^ and so is sensitive to any correction which 
shifts the thermodynamics. Therefore we particularly want an algorithm which respects 
the correct thermodynamics to good, controllable accuracy. Our choice is to implement 
Eq. (13. 7p as a series of short Hamiltonian trajectories of length t^, with the electrical field 
drawn randomly from a Gaussian ensemble at the beginning of each trajectory. To see how 
this reproduces Eq. (13. 7p . consider the following Hamilton's equations (in Aq = gauge, 
which we choose for explanatory convenience): 



dH 



OA' 
dtAi = Ei, 

{Ei{x, t = 0)Ej{y, t = 0)) = T6ij6'{x - y) . (3.8) 



(The lattice implementation of these equations is known, see for instance [18|] but note that 
in the above we do not enforce Gauss' Law.) To second order in small t^, the change in Ai 
is found to be 

A,{x, tn) - A,{x, 0) = E,{x)t^ - + 0{t^^) . (3.9) 

This is the same as Eq. (13. 7p if we interpret = t\/2. Hence, a randomization of the E fields 
followed by an evolution of time t^ corresponds to an evolution of Tl = t\/2 units of Langevin 
dynamics up to higher order in t^ corrections. Higher order corrections are inevitable in a 
discretization of Langevin dynamics, and in this case we know that any corrections to the 
interpretation of the algorithm in terms of Langevin dynamics must enter at even powers 
in since the algorithm described is actually t-even. In practice we measure Fl for two or 
more values of tn and make an extrapolation to the small tn limit, based on values no longer 
than one lattice unit. We also apply the known 0{a) lattice-continuum corrections to the 



lattice spacing 43| and to the time scale t, which is known for Hamiltonian dynamics 27 . 
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FIG. 2: (Color online) Volume dependence of the sphaleron rate in Langevin dynamics, determined 
on lattices with spacing a = 1/Ncg'^T, for SU(2) and SU(3). Errors are statistical only; data are 
not extrapolated to t = or a = 0. Using the indicated scaling, the SU(2), SU(3) behaviors look 
very similar; the rate saturates around L = Ih/NcQ^T . 

IV. WEAK-COUPLING RESULTS 



We want the large volume, small lattice-spacing limit and our goal is a few percent 
precision. There are many systematic limits to consider; lattice spacing a — 0, volume 
\^ — )■ oo, refresh frequency — ^ 0, leapfrog timestep size, and parameters of our N^^ 
measuring method. Some of these are simple to control; we can check systematically that 
our iV(^g measurement is topological because our method is "calibrated" with integrations to 



the vacuum [42|. Errors associated with the leapfrog step size are quadratic in the step; we 
compared step sizes of 0.2 and 0.1 and found a 4 ± 2% change in Fl. Hence the difference 
between 0.1 and is 1.3 ± 0.7% which is acceptably small. We typically extrapolate over 
two tn values, t^i = a and = a/2. The resulting Fl generally differ by less than 10%, so 
higher orders in the extrapolation should be at the 1% level and therefore negligible. 

Next we consider the extrapolation to infinite volume. The good news is that 3-D SU(A'c) 
gauge theory has a mass gap, so the approach to large volume behavior is expected to improve 
exponentially in box length for boxes much larger than the correlation length. In practice 
we measure the sphaleron rate at fixed lattice spacing a (technically, at a fixed value of the 
dimensionless ratio g'^aT) for several volumes to see what volume is sufficient to achieve the 
large volume limit. This was first done by Ambj0rn and Krasnitz j2o|. Our results, shown in 
Fig. [21 indicate that the volume dependence of the sphaleron rate in SU(2) and SU(3) gauge 
theory are almost identical if we scale the box length by NcQ^T, not g^T. We also see that 
the sphaleron rate increases with volume but saturates around L = Ib/N^g'^T. Therefore 
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squared latt. spacing (aN^g^T)^ squared latt. spacing (aN^g^T)^ 

FIG. 3: (color online) Finite lattice spacing extrapolation of SU(2) [left] and SU(3) [right] lattice 
Langevin sphaleron rate. The two largest lattice-spacing points were not used in the extrapolation. 



we should make sure L > Ib/N^g^T to be in the infinite volume limit. We therefore use 
L > 18 /N^g^T in what follows.^ 

Finally, we consider the small lattice spacing a limit. Since the numerical effort required 
to get the same level of statistical significance increases as a~^, we should try to understand 
this limit as well as possible. We know that, at the level of the thermodynamics of the 
system, there are corrections already at the 0{a) level. These are well understood and 



amount to an effective rescaling of the lattice spacing [43 



0.07918 SU(2) 
N^g^T X <j 0.10232 SU(3) (4.1) 

flcorrected ^bare [ 0.1208A - 1 / {6N^) SV{Nc) 

Choosing a < l/N^g'^T is then sufficient to push this correction under 10%. There are 
unknown 0{a^) corrections which are presumably of order the square of this size, so a < 
l/N^g'^T would leave of order 1% unknown corrections to the effective lattice spacing. But 
the scaling between lattice and continuum Langevin diffusion rates involves a^, so even a 
1% error in a leads to a 5% error in the diffusion rate, which is not acceptable. There are 
also other unknown 0{a'^) corrections which cannot necessarily be represented as a rescaling 
of a or of Tl. Hence we should try to get to lattices which are somewhat smaller than 
a = l/N^g'^T, and to perform an a — )■ extrapolation over the remaining corrections. 

Our data for Fl as a function of is presented in Fig. [3] and Table [H The two largest- 
a lattices used are not included in the fit, because the corrections are too large (the fit 
becomes poor) and because the behavior of A^^^^g is not yet strictly topological (the N^^ 



^ Note that using very large volumes is not helpful; since the measurement of N^^ is global over the lattice, 
there is no gain in statistics as the lattice volume increases, while the numerical effort increases with 
volume. 
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0.80 
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SU(3) 
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182.4 ± 6.0 


SU(2) 


0.67 


20.40 ± 0.70 


SU(3) 


0.75 


185.2 ±6.8 


SU(2) 


0.50 


21.56 ±0.58 


SU(3) 


0.60 


193.6 ±7.6 









TABLE I: Numerical results for the sphaleron rate under Langevin diffusion, for SU(2) and SU(3) 
gauge theory and several lattice spacings. Each result is taken at a volume with L > 18/Ncg'^T 
and is already extrapolated to zero tn] systematic errors from these extrapolations are smaller than 
the indicated statistical errors. 



Colors iVc 


Coeff. in Eq. (fOl) 


Coeff. in Eq. K2\f 


2 


9.53 ±0.14 


0.198 ± .003 


3 


127.8 ±4.0 


0.197 ± .006 


5 


2790 ± 70 


0.186 ± .005 


8 


51300 ±3100 


0.199 ± .012 



TABLE II: The Langevin equation N^^ diffusion rate as a function of the number of colors A'c, for 
a fixed lattice spacing a = 0.8/Ncg'^T and field refresh rate = a. 



measuring algorithm sometimes shows an inconsistency in which J FF from vacuum, along 
the Langevin history, and back to vacuum is far from an integer). To convert the extrapolated 
value into the result in Eq. (11. 3p . we multiply by (3iVc/47r) as indicated in Eq. (13. 3p . Eq. (13. 6p . 
As a check on our Langevin dynamics method, we have also repeated the calculation of the 
SU(2) sphaleron rate originally performed in j27|; we find that the coefficient in Eq. (II. 3p 
for SU(2) is 10.0 ± 0.3, in agreement with 10.8 ± 0.7 found there. 

In j29[ one of us speculated that Fl should scale with N^. as ^^"^^ (NcaYT^, and in 
particular 

Tsphai = ^' (^^^) (in Y + 3-0410) (N^ - l)iV,^«fT^ , (4.2) 

with K,' ~ 0.24. We have investigated this speculation by computing the sphaleron rate at a 
single value of a = 0.8/Ncg'^T for several values of A'^c, shown in Table [Tll The data in the 
table are also not extrapolated to tn = but are computed at tn = cl- Extrapolating by 
assuming that the a and scaling are the same as for Ac = 2, 3, the correction lowers 
the rate by 1.8% and the ^ correction raises the rate 9% so we find n' = 0.21 ± .01, as 
quoted in Eq. (11.61) . 
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V. FINITE-COUPLING DESPERATION 



Our calculation has made a strict weak-coupling expansion, which is really two approx- 
imations. First, in writing the theory in terms of classical fields on the lattice, we assume 
the nonperturbative IR dynamics is essentially classical. Second, by implementing Langevin 
dynamics, we assume that color conductivity is large and wavelength independent. We can- 
not relax the first approximation without losing our computational method altogether. But 
we can relax the second approximation and allow the conductivity to be finite, by studying 



the theory with Hamiltonian dynamics [18|, |20[. This allows us access to "intermediate" 
couplings, though with larger theoretical errors. 

At weak coupling, the continuum theory has three parametric length scales: 

• The length scale 1/T, where most thermal degrees of freedom reside, 

• The length scale l/y/N^gT associated with screening phenomena, 

• The length scale l/N^g'^T, where physics is nonperturbative. 

The corresponding scales for the Hamiltonian lattice theory are a, ^Ja/N^g'^T, and l/N^g'^T. 
Only two of these scales are relevant for the dynamics of topology change; the l/N^g'^T 
scale where the nonperturbative physics occurs, and the l/y/N^gT "screening" scale which 
infiuences the dynamics of the IR scale. On the lattice one does not choose a, g'^, and T 
separately; one controls the dimensionless ratio aN^g'^T (often written 2N^/f3^). But this 
choice can be paraphrased as a choice for a finite value of g'^N^, as we now show. 

Since the sphaleron rate involves the evolution of infrared magnetic fields, the relevant 



feature of 1/gT physics is the "damping rate" for the evolution of magnetic fields 3l|- In 
the continuum this is proportional to m^. On the lattice it is proportional to an integral 



which diverges as 1/a. Equating them, we find [31 



{2N, + N,)ng^T' 2.U988iTN,g^T 

«;7contin = 77Z = kjl^tt = 5 , (5.1) 

which relates the screening scales of the lattice and continuum theories. Re-organizing, we 
can determine g^ in terms of (aNcg'^T): 

2-14988 X 3 2 2^ ^ 9' 2N, + N, , , 

, AT — J^cQ = aJ^cQ J =^ — = rloAcO I ) . (5.2) 

tt^MYl 4n 12iV2(^2.14988)^ ^ ' ^ ' 

We use this relation to convert a value of the lattice spacing - really, of aN^g'^T - into an 
equivalent value of = g'^/4:7i. 



In studying the SU(2) sphaleron rate, where the gauge coupUng is smaU, it is convenient to add "extra" 



degrees of freedom to the lattice, as done in 2^l26[. This effectively shortens the screening length scale 



on the lattice, which is the same as effectively making the coupling smaller. Since we want to explore 
the largest possible couplings, this is the opposite of what we want, which is why we consider the pure 
classical lattice theory without additional degrees of freedom. 
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0.062 


0.093 


17.28 ± 0.30 


27.7 


2.00 


0.052 


0.078 


14.06 ± 0.20 


23.5 


1.72 


0.044 


0.066 


12.60 ± 0.27 


20.5 


i.OU 
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U.UOo 


ii.4o it U./u 


ICO 

io.z 


1.20 


0.031 


0.047 


10.41 ±0.22 


14.9 


1.00 


0.026 


0.039 


9.18 ±0.23 


12.7 


0.75 


0.019 


0.029 


8.51 ±0.16 


9.8 


0.60 


0.016 


0.023 


7.81 ± 0.20 


8.0 



TABLE III: SU(3) sphaleron rate under Hamiltonian dynamics, and its interpretation in terms of 
Os for the 0-flavor and 3-flavor theories. Errors shown are statistical, and are dwarfed by much 
larger theoretical errors in equating lattice and continuum theories, which we have not estimated. 
We also show the result implied by Eq. (jl.3p . which works well at the smallest lattice spacings but 
is a poor description for the coarsest lattices. 

Table IIIII presents our results for the sphaleron rate under Hamiltonian dynamics. The 
largest value of aN^g'^T corresponds to the point where we can barely distinguish topological 
behavior on the lattice. For this value, the infinite volume limit was achieved with a lattice 
8 on a side! For as larger than this corresponding value, our ability to guess the strong 
sphaleron rate gets even weaker; we must extrapolate based on the behavior shown in the 
table. This corresponds to the case relevant for heavy ion collisions. The figure also shows 
the performance of the leading-log expansion, that is, the guess for the sphaleron rate using 
Eq. (11. 3p together with the matching for the Debye screening scale from Eq. (15. 2p . The 
leading-log (Bodeker effective theory) approach reproduces the lattice Hamiltonian result 
for the finest lattices but is rather far off for the coarsest lattices, showing that the evolution 
is not overdamped for such coarse lattice spacing (large ttg). 

VI. DISCUSSION 

We have computed the sphaleron rate in SU(3) gauge theory at weak coupling. In the 
weak coupling limit, there is a rigorous relation between the sphaleron rate and the topo- 
logical diffusion rate for classical lattice gauge theory under Langevin dynamics, which we 
have computed with few-% error bars. We have also shown that the A^^c-dependence for the 
sphaleron rate has a surprisingly simple and intuitive behavior. Our analysis of Hamiltonian 
dynamics on coarser lattices suggests that this Bodeker effective description only works at 
very small values of a^- While it should apply for the SU(2) electroweak sector and for QCD 
at GUT scale temperatures, it is not a successful description at electroweak temperatures. 

For applications to electroweak baryogenesis, we need the strong sphaleron rate in an 
A^f = 6 theory with as ~ 0.1. The value of as for A'f = 6 is twice the value for A'f = 0, so 
this corresponds roughly to our [aN^g'^T) = 2 data point in Table IIII( the strong sphaleron 
rate relevant at the electroweak scale is approximately 

rstrongsphal(iVf = 6, C^s ~ 0.1) ~ 14a^T^ . (6.1) 
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It is difficult to assign theoretical error bars to this estimate but they should be rather large; 
the relevant lattice spacing is so large that the topological nature of the FF measurement 
is barely under control, and the lattice spacing corrections discussed earlier are not small. 
It is also not clear that the lattice field dynamics are a good description of the continuum, 
not-so-small coupling dynamics. We only have rigorous arguments that they agree for small 
coupling.^ The continuum dynamics have 0{g) (technically 0{Ncg'^T/mj^)) corrections, the 
lattice dynamics have 0(\/ Nr.q'^aT) corrections, and there is no guarantee that they are the 
same (see for instance [4^). Therefore we expect at least (9(50%) systematic corrections to 
this result. 

Unfortunately, the coupling values ttg ~ 1/3 relevant for heavy ion physics correspond to 
lattices too coarse for any study of topology to be successful. Therefore the best we can do 
is to extrapolate our existing results towards this regime. One might guess based on Table 
Hm that the sphaleron rate at large couplings is Fsphai ~ SOa^T^. But this is a desperate 
extrapolation. It is not clear to us how to proceed at such large couplings; as we have 
argued, analytic continuation from Euclidean lattice calculations seems to be even more 
difficult than for other transport coefficients. But at least the behavior at the electroweak 
temperature scale, and at higher temperatures, is under some control. 
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